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Any injection of electromagnetically interacting particles during the cosmic dark ages will lead to 
increased ionization, heating, production of Lyman-a photons and distortions to the energy spec¬ 
trum of the cosmic microwave background, with potentially observable consequences. In this note 
we describe numerical results for the low-energy electrons and photons produced by the cooling 
of particles injected at energies from keV to multi-TeV scales, at arbitrary injection redshifts (but 
focusing on the post-recombination epoch). We use these data, combined with existing calculations 
modeling the cooling of these low-energy particles, to estimate the resulting contributions to ioniza¬ 
tion, excitation and heating of the gas, and production of low-energy photons below the threshold for 
excitation and ionization. We compute corrected deposition-efficiency curves for annihilating dark 
matter, and demonstrate how to compute equivalent curves for arbitrary energy-injection histories. 
These calculations provide the necessary inputs for the limits on dark matter annihilation presented 
in the accompanying Paper I, but also have potential applications in the context of dark matter de¬ 
cay or de-excitation, decay of other metastable species, or similar energy injections from new physics. 
We make our full results publicly available at http://nebel ,rc .fas.harvard.edu/epsilon, to facil¬ 
itate further independent studies. In particular, we provide the full low-energy electron and photon 
spectra, to allow matching onto more detailed codes that describe the cooling of such particles at 
low energies. 

PACS numbers: 95.35.+d,98.80.Es 


I. INTRODUCTION 

Between recombination and reionization, the universe 
experienced an epoch of extremely low ionization, known 
as the “cosmic dark ages”. If new physics were to in¬ 
ject electromagnetically interacting particles into the uni¬ 
verse during this period - with the classic examples being 
dark matter (DM) annihilation or decay - and conse¬ 
quently induce increased ionization, it could broaden the 
last scattering surface and have striking effects on the 
anisotropies of the cosmic microwave background (CMB) 
m- Furthermore, heating of the gas induced by such 
energy injections could have observable effects on the 
21cm line from neutral hydrogen 0HZ], and the produc¬ 
tion of additional low-energy photons could distort the 
blackbody spectrum of the CMB MU- 

A critical question in studies of the observational con¬ 
sequences of energy injection is what fraction of the in¬ 
jected power proceeds into the various observable chan¬ 
nels, and over what period of time. As discussed in 
S3, photons and e + e pairs injected around the elec- 
troweak scale (a typical scenario in annihilating DM 
models) promptly convert the bulk of their energy into 
photons with energies lying within a redshift-dependent 
semi-transparent “window” [2], where the dominant cool¬ 
ing mechanisms have timescales comparable to a Hubble 
time. Some fraction of these photons never scatter again, 
slowly redshifting and contributing to X-ray and gamma- 
ray background radiation in the present day; others even- 
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tually partition their energy into lower-energy photons 
and electrons, which are either efficiently absorbed by 
the gas or contribute to distortion of the CMB spectrum. 
Accordingly, energy may be deposited and contribute to 
observable signatures at times long after its original in¬ 
jection. 

In this work we employ the code initially described in 
m and refined in mm to describe the energy depo¬ 
sition histories corresponding to particle injection at ar¬ 
bitrary energies and redshifts. In m we computed the 
partition between “deposited” energy and free-streaming 
high-energy photons, for injections of photons, electrons 
and positrons at arbitrary energy and redshift. In that 
work “deposited” energy was taken (as in J12]) to en¬ 
compass low-energy particles in general, including dis¬ 
tortions to the CMB spectrum and ionization, excitation 
and heating of the gas. 

In order to convert from this overall “deposited en¬ 
ergy” to the individual deposition channels, articles in 
the literature have generally followed [2] in employing a 
simple prescription for the fraction of deposited power 
proceeding into ionization, excitation and heating, based 
on studies of this fraction for 3 keV electrons m- More 
careful modeling of the cooling of electrons and photons 
had supported this estimate mmm, in particular for 
the fraction of deposited power proceeding to ionization 
(which is the most important channel for determining 
the impact on the CMB anisotropy spectrum, e.g. W- 
However, m demonstrated that this prescription can be 
quite inaccurate in general, as somewhat higher-energy 
electrons (between a few keV and a few MeV in kinetic 
energy) deposit the bulk of their energy into distortions 
of the CMB spectrum rather than through interactions 
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with the gas. By employing the code developed for [T2] 
to model the cooling of high-energy particles (and the 
secondary particles produced by their cooling) down to 
3 keV energies, and then matching onto a separate code 
handling the cooling of electrons below 3 keV, m and 
[20] presented updated estimates of the power proceed¬ 
ing into ionization - and hence the constraints from the 
CMB - for a selection of DM models. 


In this note, we present a similar update for all injec¬ 
tions of photons and e + e _ pairs at redshifts during the 
cosmic dark ages, with injection energies in the 0(keV- 
TeV) range. 1 Furthermore, we provide estimates for the 
power proceeding into Lyman-a photons, heating, and 
CMB spectral distortion (by continuum photons below 
10.2 eV), as well as ionization. We provide the full spec¬ 
tra of photons and electrons below 3 keV produced by the 
high-energy code at all timesteps, for all injection ener¬ 
gies and redshifts, to facilitate the interfacing of these 
results with more detailed and precise models for the 
cooling of the low-energy particles. 

In Section in! we review the issues that mandate an 
improved treatment of the energy deposition, with a fo¬ 
cus on setting constraints on energy injection via the 
anisotropies of the CMB. In Section [TTT| we review the key 
elements of the code employed, and describe the resulting 
publicly available dataset. In Section IV we review two 
procedures for converting the low-energy spectra into es¬ 
timates for the deposition to various channels, and show 
updated results for the total energy deposited into the 
various channels under these prescriptions. In Section 
|Vj we review how to determine so-called f(z ) curves - 
the power deposited at any given redsliift, normalized to 
the power injected at that same redshift - for any en¬ 
ergy injection history; as an example, we present f(z) 
curves corrected for the systematic effects identified in 
El. for general DM annihilation models, suitable for 
use with studies that employed earlier injection-energy- 
independent prescriptions for the fraction of deposited 
power proceeding into ionization. Finally, we present our 
conclusions. Appendix [A] provides detailed descriptions 
for the files containing our results, available online at 
http://nebel.rc.fas.harvard.edu/epsilon. 


II. CONTINUUM PHOTON LOSSES AND 
ENERGY DEPOSITED TO THE GAS 


To determine the constraints on any model of new 
physics that injects electromagnetically interacting parti¬ 
cles into the universe during the cosmic dark ages, the key 
figure of merit is the power deposited into the relevant 


channel(s) at any given redshift. For example, for con¬ 
straints based on the anisotropies of the CMB, the most 
important channel is ionization, 2 and so the constraints 
are determined by the power deposited into ionization 
of the gas as a function of redshift. The distortions to 
the CMB energy spectrum and the gas temperature are 
non-zero, but the constraints arising from those channels 
are much weaker than from the impact of extra ioniza¬ 
tion on the CMB anisotropies (e.g. 0). Once computed, 
the power deposited as extra ionization can be incorpo¬ 
rated into public codes describing recombination - such 
as RECFAST [23j, CosmoRec (23] or HyRec [25] - as de¬ 
scribed in mm, and the resulting ionization history used 
to compute the effects on the CMB. 

Prior studies (e.g. m 0 cca EDi EM]) have divided 
the calculation of the power deposited to each channel 
into two steps: (1) computing the total deposited power 
as a function of redshift, and (2) computing the fraction 
of deposited power that proceeds into each channel. The 
result of step 2 has been presumed to be a function of the 
gas ionization fraction only, independent of the details 
of the energy injection; under this assumption, all the 
dependence on the energy injection model is partitioned 
into step 1. The results of step 1 have frequently been 
approximated as a constant efficiency factor /, so that 
the deposited power at any redshift is simply /x the 
power injected at that redshift, and / captures all model- 
dependence in how the energy is deposited (e.g. [51 [27]). 

This approach assumes that the fraction of deposited 
power proceeding into each channel is independent of the 
spectrum of particles marked “deposited”; that once the 
energy contained in free-streaming high-energy photons 
has been removed, only the total energy of the remaining 
particles matters. As shown in El. this approach can 
lead to constraints that are incorrect at the factor-of-two 
level (even within the limited parameter space explored 
there). 

The power into ionization has frequently been esti¬ 
mated in the literature using a simple ansatz first devel¬ 
oped by Chen & Kamionkowski 0, based on an earlier 
numerical result for ionization by 3 keV electrons m ■ 
(1 — x e )/3 of the absorbed power goes into ionization, 
where x e is the ambient hydrogen ionization fraction. 
We denote this ansatz as “SSCK” (Shull, van Steenberg, 
Chen & Kamionkowski). This estimate has been sup¬ 
ported by more recent and detailed calculations m na¬ 
ns]. M studied the effect of taking into account the spec¬ 
trum of electrons below 3 keV, as opposed to simply using 
the results for 3 keV electrons, in determining the fate of 
the deposited energy; the effect was found to be small for 


1 We do not in this work provide a detailed study of the energy 

losses of protons and antiprotons; an approximate method for 
including these contributions to energy deposition can be found 
in [2T1 . 


2 There is a subdominant effect from the production of additional 
Lyman-a photons, since atoms in an excited state can be more 
easily ionized by the ambient CMB photons; however, neglecting 
this effect entirely has been shown to change the constraints at 
only the 5% level EHEal, justifying a simplified approximate 
treatment of this contribution. 
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the CMB constraints on the DM models studied in that 
work. The reason is that the fraction of power deposited 
as ionization is relatively stable for ~ 100 eV - 3 keV 
electrons; we show this fraction as a function of redshift 
(employing a background ionization history from REC- 
FAST as in [34| . with no new energy injection) in Figure 
|TJ for a range of electron energies. We also show the esti¬ 
mate from the “SSCK” energy-independent prescription. 
We will refer to the approach of taking the results pre¬ 
sented in m for 3 keV electrons, and applying them to 
estimate the fraction of deposited power proceeding into 
the different channels, as the “3 keV” prescription; this 
is the approach employed by the Planck Collaboration in 
setting constraints on annihilating DM [33]. 



i+z 


FIG. 1: Fraction of deposited power that contributes to (hy¬ 
drogen) ionization as a function of redshift, employing the 
“SSCK” prescription (solid black line), and the results pre¬ 
sented in m for electrons at a range of energies. The “3 
keV” prescription corresponds to the purple dashed line. 

However, as noted in [14j . higher-energy electrons (at 
kinetic energies of a few keV to a few MeV) lose a very 
large fraction of their initial kinetic energy to inverse 
Compton scattering on the CMB; for mildly relativis¬ 
tic electrons, the resulting upscattered photons have too 
little energy to further interact with the gas after recom¬ 
bination (being well below the excitation and ionization 
thresholds). Accordingly, the amount of energy going 
into ionization, excitation or heating of the gas is sup¬ 
pressed, relative to the case for 3 keV electrons. This 
effect was underestimated in some earlier studies USEE] 
due to a mistake in the expression for the cooling time 
due to inverse Compton scattering, for non-relativistic 
electrons (the correct expression is given in e.g. T2|). 
Consequently, for particles injected at high energy, the 
details of the low-energy electron and photon spectra 
produced by their cooling can significantly influence the 
fraction of deposited power that proceeds into ionization. 

Equivalently, describing power degraded to low energy 
scales as “deposited” can be somewhat misleading, since 
photons at energies comparable to the CMB are not ab¬ 
sorbed by the gas. The two-step approach described 
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Electron energy (eV) 

FIG. 2: Fraction of electron kinetic energy eventually par¬ 
titioned into (red/pink) excitation/ionization/heating of the 
gas, plus electrons below the energy cutoff where atomic cool¬ 
ing processes dominate, and (blue/violet) photons produced 
by inverse Compton scattering of the electron on the CMB. 
Photons below 10.2 eV, which no longer interact significantly 
with the gas, constitute a subset of the latter contribution 
(light/dark green). All curves are for electrons depositing 
their energy at redshift 600. Dashed and dotted lines (pink, 
blue, light green) correspond to a 1 keV threshold, whereas 
solid lines (red, violet, dark green) correspond to a 3 keV 
threshold; we see that the behavior is fairly independent of 
the chosen threshold except very close to it. 

above could be improved by redefining “deposited” en¬ 
ergy to exclude photons below 10.2 eV as well as free- 
streaming high-energy photons. However, such photons 
are also produced by the cooling of ~ keV and lower- 
energy electrons; thus with this definition there is no 
range of energies within which particles can be treated 
as contributing solely to “deposited energy”. 

Figure [2] demonstrates the magnitude of the fractional 
energy loss to photons too low-energy to excite hydrogen, 
as a function of the initial electron energy, at a redshift of 
z; = 600 (where DM annihilation typically has its greatest 
effect on the CMB anisotropies |M|). Electrons that cool 
to some threshold (set to 1 keV or 3 keV) are presumed 
to lose all their energy to atomic processes, which dom¬ 
inate inverse Compton scattering in that energy range 
and are well-described by existing low-energy codes (for 
example, the fraction of power proceeding into ioniza¬ 
tion is well-characterized by the curves shown in Figure 
[l]). An electron above this threshold loses some fraction 
of its energy to atomic processes and some to inverse 
Compton scattering, in the process of cooling down to 
the threshold, producing photons on the way. As shown 
in the figure, when electrons are injected with (kinetic) 
energies between a few keV and ~ 10 MeV, the bulk of 
its energy is lost into photons below 10.2 eV in energy, 
for which the universe is approximately transparent. 

It is difficult to properly model these inverse Compton 
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losses in the context of Monte Carlo simulations for the 
low-energy atomic processes, due to the huge number of 
nearly-elastic collisions that are involved. As discussed 
above, this difficulty was avoided in m by using two sep¬ 
arate codes to treat the low-energy (below 3 keV) and 
high-energy (above 3 keV) electrons. The high-energy 
code degrades the initial particles from their injected en¬ 
ergy down to the 3 keV threshold, fully taking into ac¬ 
count the effects of both inverse Compton scattering and 
redshifting. The results of the high-energy code - spectra 
of electrons and photons below 3 keV - can then be fed 
as inputs into a low-energy Monte Carlo code that treats 
the complex atomic cooling processes in detail. We take 
the same approach here, and make the results public. 


III. DESCRIPTION OF NUMERICAL RESULTS 
A. Review of the numerical method 

The code developed for |12| takes as an input some in¬ 
jection of photons and electrons, with a specified redshift 
and energy dependence. Backreaction on the CMB pho¬ 
tons and gas is not included, as large modifications to 
the ionization history or CMB spectrum are ruled out by 
observational constraints, and consequently the problem 
is (to a good approximation) linear. We thus populate 
individual energy bins with electrons/positrons or pho¬ 
tons at a specific redshift and track the spectral evolution 
with redshift. Our 40 energy bins are log-spaced between 
1 keV and 10 TeV, in photon energy and electron kinetic 
energy. We employ 63 log-spaced redshift bins spanning 
the range from 1 + z = 10 to 3000 (for the redshift at 
which the energy is injected - for the redshift at which 
the energy is deposited, more finely-binned results are 
available upon request). 

At each timestep, the photon spectrum is updated with 
the results of the various scattering and pair produc¬ 
tion processes described in [12] . and redshifted. Photons 
at sufficiently low energies are tagged as “deposited”, 
stored, and removed from the part of the code that de¬ 
scribes redshifting. The threshold for this “deposition” 
occurs when the photon would on average photoionize 
an atom once per timestep, and as in m, we choose a 
timestep of dln(l + z) = 1CU 3 (it was confirmed in that 
work that the results were converged at such a timestep). 

As described in [12] , the free electrons produced at each 
timestep (by direct injection, pair production cascades, 
Compton scattering, etc) lose their energy on timescales 
much shorter than a Hubble time, and are handled by 
a separate module that includes inverse Compton scat¬ 
tering and atomic cooling processes, and resolves their 
energy deposition entirely at each timestep. 

We choose a threshold of 3 keV as the separation scale 
between “low” and “high” energies. At “low” energies, as 
described above, inverse Compton scattering of electrons 
and redshifting of photons can be neglected and the re¬ 
sults of a detailed Monte Carlo code will be employed to 


model electron cooling through atomic processes. “De¬ 
posited” photons above 3 keV (that is, efficient photoion¬ 
izers that will produce electrons above the threshold) are 
assumed to immediately photoionize hydrogen, convert¬ 
ing their energy into a spectrum of secondary electrons 
(the sub-percent energy loss to the ionization itself is ne¬ 
glected for these photons). “Deposited” photons below 3 
keV are not processed further by the high-energy code, 
but are stored as an output at every timestep, to facilitate 
their use as an input for more detailed codes describing 
the low-energy cooling. Photons below 3 keV are only 
entirely tagged as “deposited” if the timescale for pho¬ 
toionization is less than 10 _3 x a Hubble time (for that 
photon energy and redshift); otherwise, only the frac¬ 
tion of the photons that would photoionize the gas in 
that timestep are tagged as “deposited” and added to 
the output (since their photoionization will create sec¬ 
ondary electrons below 3 keV), while the remainder are 
tracked in the main code. 

Electrons above the 3 keV threshold lose energy domi¬ 
nantly through inverse Compton scattering, as discussed 
above, but simple estimates for energy losses to exci¬ 
tation, ionization, and heating by electrons above the 
threshold are stored as outputs to the code (and consti¬ 
tute a small contribution to the overall excitation, ioniza¬ 
tion and heating due to the choice of threshold, typically 
~ 5 — 15 % or less). These estimates are based on the 
implementation of the excitation, ionization and heating 
cross sections described in |12j . which have been con¬ 
firmed to give results fairly consistent with the detailed 
low-energy code m , with differences at the level of a few 
percent. Calculations which treat the low-energy cooling 
in detail can separate Lyman-a photons from sub-10.2 
eV photons produced by collisional excitation, but for 
the subdominant contribution from the high-energy code, 
we simply assign all the power deposited via excitation 
to the Lyman-a channel. 

Sub-threshold photons produced by inverse Compton 
scattering of these above-threshold electrons are also 
stored as outputs; since for mildly relativistic electrons 
the boost from inverse Compton scattering is small, we 
must also take into account the removal of the original 
photons from the CMB. (That is, if a 1 eV photon is up- 
scattered to 2 eV, it is described as the removal of 1 eV of 
energy from the CMB and the addition of a 2 eV photon 
to the output photon spectrum.) For the electron ener¬ 
gies where the impact on the CMB is non-negligible, the 
inverse Compton scattering cross section is independent 
of the initial photon energy, so the removal of photons 
does not change the spectral shape and it is sufficient to 
characterize it in terms of energy loss. Once electrons 
cool to the 3 keV threshold or below (electrons produced 
by photoionization, Compton scattering or pair produc¬ 
tion may be produced with less than 3 keV of kinetic 
energy), they are added to the output low-energy elec¬ 
tron spectrum for that timestep. 

As in m , we are primarily interested in charge-neutral 
sources of energy injection, so we will generally consider 
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injecting e“*"e _ pairs rather than electrons alone. This is 
largely irrelevant for the energy-loss mechanisms impor¬ 
tant for high-energy particles, but when the positrons 
have cooled far enough, they annihilate with ambient 
electrons, producing gamma rays. For a relativistic 
positron, the vast majority of its energy will be deposited 
via inverse-Compton-scattered photons and their sub¬ 
sequent cooling, and these processes do not depend on 
the charge of the original particle; for a non-relativistic 
positron, the bulk of its energy will be contained in its 
mass energy, and be deposited via the photons from an¬ 
nihilation. Accordingly, aside from tracking the annihila¬ 
tion photons we do not distinguish between positrons and 
electrons; since the low-energy electrons produced domi¬ 
nantly arise from photoionization and Compton scatter¬ 
ing, they are indeed electrons (not positrons) and can be 
treated as such for detailed low-energy codes. Injection 
of electrons without accompanying positrons can there¬ 
fore be modeled using the results for e + e“ pairs, together 
with the results for photon injection to remove the im¬ 
pact of the photons from annihilation (see m for an 
example of the procedure). 

We truncate the calculation at z = 9, and advise cau¬ 
tion in using our results for z < 30, due to two simplifying 
assumptions in the code 3 : first, no model for reioniza¬ 
tion is included in our baseline ionization history, and 
second, interstellar radiation fields other than the CMB 
are not included, and may be relevant after the onset of 
star formation. To properly include these effects would 
involve considering a range of models for ionization and 
star formation, and testing the sensitivity of the results 
to these choices; we thus defer a detailed study of low 
redshifts to future work. The effect of reionization on 
the cooling of high-energy electrons and photons is likely 
to be rather small, as inverse Compton scattering and 
the main energy-loss processes for high-energy photons 
are largely insensitive to the ionization fraction |T2t, so 
one might expect the output spectra of low-energy pho¬ 
tons and electrons described in this section to be fairly 
independent of the choice of reionization model. How¬ 
ever, the rate at which low-energy photons are absorbed 
by photoionization, and the partition of the low-energy 
photons’ and electrons’ energy in to t he various deposi¬ 
tion channels (described in Section IV), will both depend 
strongly on the background ionization fraction. 


B. Structure of the outputs 

The results of the high-energy code are expressed as 
entries in a three-dimensional grid for each of the particle 
types, describing - for a particular redshift-of-injection, 
initial energy (for e + e - pairs, this corresponds to the 
initial kinetic energy of one member of the pair), and 


3 We thank Aaron Vincent for raising this question. 


redshift-of-deposition - the following outputs: 

• An estimate of the energy deposited by electrons 
above 3 keV in this timestep, as a fraction of the 
injected energy (for e + e _ pairs, this includes the 
mass energy), into the following channels: 

1. Ionization of the H gas, 

2. Ionization of the He gas (set to zero for this 
high-energy deposition, but listed for com¬ 
pleteness as this channel will be populated 
later), 

3. Excitation of the gas / production of Lyman-a 
photons which can excite neutral hydrogen, 

4. Heating of the gas, 

5. Production of photons with insufficient energy 
to either excite or ionize the gas (i.e. distor¬ 
tion of the CMB spectrum). 

• The total energy of photons removed from the CMB 
by upscattering, in this timestep, as a fraction of 
the injected energy. 

• The spectrum of low-energy electrons (below 3 
keV) produced in this timestep (and the corre¬ 
sponding array of energy bins), expressed as the 
spectrum dN/dE of electrons per pair of injected 
particles (i.e. in the case where DM annihilation is 
the source of energy injection, this is equivalent to 
the spectrum per annihilation). 

• The spectrum of low-energy photons (below 3 keV) 
produced in this timestep (and the corresponding 
array of energy bins), expressed as the spectrum 
dN/dE of photons per pair of injected particles. 

As a reminder, all timesteps have dln(l + z) = 10 -3 ; 
each of these results can be divided by this quan¬ 
tity to obtain (approximately) timestep-independent 
functions for the energy deposition and low-energy 
spectra. These results are available online at 
http://nebel.rc.fas.harvard.edu/epsilon, in the 
form of FITS and .dat files, and are described in detail 
in Appendix [A] 

As outlined in Table [T] we label the resulting spectra 
and energy injections by S s ^/ es or . The label 

“species” is either “ 7 ” (for injected photons) or “e + e _ ” 
(for injected pairs). The label c runs from 1 to 5 and in¬ 
dexes the absorption channels: ionization on hydrogen, 
ionization on helium, excitation, heating, and production 
of photons too low-energy to interact with the gas. Since 
the direct energy absorption from high-energy electrons 
is generally subdominant (< 15%), we do not distinguish 
between ionization on hydrogen and helium from these 
electrons, assigning the ionization contribution entirely 
to channel c = 1. Due to numerical issues (truncation of 
the energy binning for both photons and electrons), the 
power lost from electrons due to inverse Compton scat¬ 
tering is not identical to the power gained by (tracked) 
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photons (the former quantity is generally slightly larger); 
the difference is assigned to channel c = 5, so it will later 
be added to the power stored in low-energy continuum 
photons, and verified to be small. The depletion of pho¬ 
tons from the CMB spectrum is tracked in the 
array; when we compute the partition of deposited energy 
into each channel, this array will be subtracted from the 
fifth channel 5g I )®^ les . 

The i,j, k labels index the redshift of deposition, en¬ 
ergy of injection and redshift of injection respectively. 
The “sec” label indicates the species of low-energy (below 
3 keV) secondary particles being described, and l indexes 
the energy of these secondaries (or kinetic energy, in the 
case of electrons). 

These data hold all the key information from the high- 
energy code, and we recommend their use as inputs for 
detailed studies of the electron and photon cooling at 
low energies. However, for ease of use we also provide 
estimates of these quantities based on the coarsely-binned 
results for the partition of electron energy given in nn. 
following the methods outlined in that work. 


IV. COMPUTING THE DEPOSITED ENERGY 
BY CHANNEL 


A. The low-energy photons and electrons 


Our next goal is to estimate how low-energy electrons 
and photons lose their energy, and thus convert the de¬ 
rived low-energy electron and photon spectra (those in¬ 
cluded in the arrays) into contributions to the 

five channels described above. These contributions can 
then be added to those obtained directly from the high- 
energy code (the latter are generally subdominant). Our 
final result will be a three-dimensional grid for 

each channel c and injected species, where as previously 
i,j and k respectively index the redshift of deposition, en¬ 
ergy of injection and redshift of injection. The elements 
of this array will correspond to the energy deposited to 
channel c in the relevant timestep for deposition, as a 
fraction of the total injected energy. In general we have: 


ry-ispecies _ ospecies 

1—4 ,ijk ^1—4 ,ijk 

+ contributions from S^cJjM 

^species _ ^species ^species 

"o ,ijk u 5,ijk ini 


"moss ,ijk 


+ contributions from S 


(1) 


For the purposes of this work, we will use simplified 
results from existing Monte Carlo codes that model the 
atomic cooling processes, as presented in m- However, 
these results are only directly applicable to low-energy 
electrons. For low-energy photons (below 3 keV) above 
13.6 eV, we follow the “best” method outlined in M- 


• For “deposited” photons with 13.6 eV < E < 3 
keV, we assume prompt photoionization leading to 


a secondary electron and a contribution (of 13.6 eV 
per ionization) to channel c = 1. It was shown in 
1 14] that separating out ionization on helium has a 
negligible impact on the results. This contribution 
is given by: 


^-.species 

l,ijk 


E 

l> 13.6eV 


13.6eV 


S^Ej ec d\nE sec 
2 {Ei + m s P ecies ) ’ 


( 2 ) 


where the denominator is the total injected energy 
(i.e. twice the injection energy for photons, or twice 
the (injection kinetic energy + electron mass) for 
e + e _ pairs). 

• For photons with 10.2 eV < E < 13.6 eV, we assign 
the associated energy to channel 3, since such pho¬ 
tons cannot ionize neutral hydrogen, but can excite 
it to a state from which it can be more readily ion¬ 
ized. This contribution is given by: 


^^-ispecies _ 


E 


10.2eV<Z<13.6eV 


S^(Ej ec ) 2 d]nE sl 
2 {Ei + m s P ecies ) 


(3) 


• For photons with E < 10.2 eV, we assign the asso¬ 
ciated energy to channel 5: 


species _ 


E 

Z<10.2eV 


S^;lT(EL) 2 d^E s , 

2 {Ei + TO s P ecies ) 


(4) 


This approach may be inaccurate at redshifts signifi¬ 
cantly prior to recombination, when the photoionization 
rate is highly sensitive to the small neutral fraction, and 
prompt photoionization may not be possible (or may only 
be possible on helium). However, small changes to the 
gas ionization fraction at these same redshifts do not af¬ 
fect the CMB anisotropies (since the universe is opaque 
to CMB photons); accordingly, this inaccuracy does not 
impact CMB constraints (see e.g. [20] where this was 
explicitly tested). After recombination, the photoioniza¬ 
tion rate is always fast relative to the Hubble time for 
some range of (low) photon energies [12] , 

This procedure converts the low-energy photon spec¬ 
trum to a low-energy electron spectrum, plus contribu¬ 
tions to channels 1, 3 and 5. We simultaneously define 
the related quantity, 


rpspecies ^species ^species 

^ corr ,ijk ^5 ,ijk ^loss ,ijk 


+ E 

Z<10.2eV 


S^{Ej ec )^\nE sec 
2 (Ei + m s P ecies ) 


(5) 


which describes the net energy proceeding into channel 
5 omitting the direct contribution from low-energy elec¬ 
trons (so this quantity can be used to correct the results 
of previous studies, which accounted only for the latter). 

To characterize the fate of the low-energy electrons, 
we employ the results presented in HU, based on Ha 
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Definition 

Channel c 

Hydrogen ionization: c = 1 or “ionH” 

Helium ionization: c = 2 or “ionHe” 

Excitation: c = 3 or “exc” 

Heating: c = 4 or “heat” 

CMB spectral distortion or “continuum” photons: c = 5 or “cont”. 

“ion” indicates the sum of the “ionH” and “ionHe” channels. 

“corr” indicates the contribution to channel c — 5 from the low-energy photons produced at each 
time-step by the cooling of high-energy particles (as opposed to the cooling of sub-3-keV electrons); this 
term was not properly taken into account in earlier prescriptions. 

“loss” indicates the depletion of the CMB by scattering (see text). 

S speclea {z,E,z') 

describes the dimensionless rate (normalized to injected energy, and differential with respect to 
dln(l + z)) at which energy is deposited into channel c, via the interactions of (primary and secondary) 
particles as they cool down to 3 keV (kinetic) energy, for a particle of the indicated species injected at 
redshift z' and energy E (see text for details). 

^species 

u c,ijk 

^.species(£* _£.?, ,jW)c(ln(l + «d ep ) = discretized form of the function above, for a specified grid of 

injection/absorption redshifts and injection energies {z^ ep , E\ z^}. 

S s s p r ea (z,E,z',E aec ) 

describes the rate dE d 'din(i+z) at which sub-3-keV secondary photons (electrons) of (kinetic) energy E sec 
are produced (“sec” labels the species of the secondaries), by the interactions of (primary and secondary) 
particles as they cool down to 3 keV (kinetic) energy, for a particle of the indicated species injected at 
redshift z' and energy E. 

ospecies 
*~sec ,ijkl 

S’Iec ecles (- z dep: F\ zF , El ec )d ln( 1 + z l Aep ) = discretized version of the function immediately above, for a 
specified grid of injection/absorption redshifts, injection energies and secondary-particle energies 
r pt yk pi \ 

L^depf J - / i ^inj ? -^sec / • 

T apec ' ea (z,E,z') 

describes the dimensionless rate (differential with respect to dln(l + z)) at which energy is absorbed into 
channel c, for a particle of the indicated species injected at redshift z' and energy E (see text for details). 
Derived from the S functions described above. 

rr\ species 
± c,ijk 

/■pspeoes( 2 :^ ep , E^ , zjC )d In(1 + 2de P ) = discretized form of the function immediately above, for a specified 
grid of injection/absorption redshifts and injection energies {z^ ep , E^ , z^ ep }. 

fc(z) 

(energy deposited to channel c in a redshift interval dz ) / (energy injected in the same interval dz) 

/high,c(2) 

as fc(z), but only including energy deposited by electrons as they cool down to 3 keV. 

/(*) 

Sc fc{z) = (total energy absorbed to all channels in a redshift interval dz) / (energy injected in the 
same interval dz) 

Xc(z) 

fc{z)/ f(z) = fraction of total absorbed energy at redshift z proceeding into channel c. 


fraction of total absorbed energy at redshift z proceeding into channel c, under an earlier simplified 
prescription labeled by “base”, corresponding to either the “SSCK” or “3 keV” prescriptions (see text). 

/ C ’ baSe (z) 

/c( 2 )/Xc ase ( x ) = should be used to replace f(z) in analyses where “base” prescription was assumed and the 
signal is determined by channel c. Would be equal to f(z) ( not f c {z)) if “base” prescription were correct. 

r n \z) 

f(z)( 1 — Xcorr(zO) = simplified estimate for / c,base (z), for c = 1 — 4, that is independent of channel c and 
prescription “base”. Can be used to estimate the corrected constraints/detectability for signals that 
depend on a combination of ionization, excitation and heating. 

F aec (z, E sec ) 

(spectrum dN/dE of secondary particles, with species labeled by “sec”, produced in a redshift interval dz) 

/ (pairs of primary particles injected in the same interval) 


TABLE I: Definitions of the various functions and index labels used in this article. Those entries above the double line describe 
completely general energy injections, whereas those below the double line require specification of an energy injection history. 


and with refinements as described in [18] n~9l . These 
data describe the energy deposition fractions by chan¬ 
nel, for electrons with energies from 14 eV to 10 keV, for 
background ionization fractions ranging from 10~ 4 to 1. 
These results are partially reproduced in Figure [l] show¬ 
ing the fraction of power proceeding into ionization as a 
function of redshift and electron energy. 

We interpolate these results logarithmically in both en¬ 
ergy and background ionization fraction. The energy 
deposition is partitioned into ionization on hydrogen, 
ionization on helium, Lyman-a photons, heating, and 
continuum photons. Integrating our low-energy electron 


spectrum over this transfer function at each triple of in¬ 
put redshift, injection energy and output redshift, we 
obtain the remaining contributions to all five channels 
c = 1 — 5. To evaluate the transfer function, we as¬ 
sume the baseline ionization history given by RECFAST 
(as shown in e.g. [34j). since the energy partition is not 
sensitive to changes in the ionization fraction m at the 
©(lO -4 ) level that is currently allowed by observations 
(e.g. [2Q] S3]). Specifically, if the background ionization 
fraction is x e and the (interpolated) fraction of energy 
deposited to channel c by an electron of kinetic energy 
E is F c (x e ,E), the contribution to a channel c from the 





















cooling of these low-energy electrons is given by: 

A Jispecies _ d\nE sec ( E l sec ) S ed jkiF c (x e (z l ) , E l sec ) 
c ,ijk — 2(Ei + m. species ) 

(6) 

This completes our main calculation of the T^ c k cs 
grid. 

1. An alternate simplified method 

Alternatively, and as a cross-check on our method, we 
can employ the “approx” procedure described in |l4j . In 
this approach, rather than compute the energy proceed¬ 
ing into each deposition channel separately, we note that 
the large corrections relative to the “3 keV” and “SSCK” 
prescriptions are driven by energy losses into continuum 
photons produced by the cooling of high-energy elec¬ 
trons, which are not accounted for by the modeling of 
3 keV electrons. This contribution to channel 5 is stored 
in the array T^ff ^ k , as defined in Equation [ 5 J If this 
contribution is subtracted from the overall deposited en¬ 
ergy, then the fractions of the remaining deposited en¬ 
ergy proceeding into channels 1-4 can be approximated 
by the “SSCK” or “3 keV” prescriptions. If the fraction 
of deposited energy proceeding into a particular channel 
c is Xc ase ( z ) i n these prescriptions, where “base” can be 
“SSCK” or “3 keV”, this approach leads to an alternate 
T c,ijk S rld S lven by: 



(7) 

This provides a reasonable approximation for channels 
1-4 but should not be used for channel 5 (one approach 
is to re-add the T^^ s k contribution to channel 5). 

B. Results for the deposited energy by channel 

As was previously done for the overall energy deposi¬ 
tion EE we can ask e.g. what fraction of a particle’s 
energy is eventually deposited to each of these five chan¬ 
nels; alternatively, upon specifying a redshift history for 
the energy injection we can ask how much power is de¬ 
posited to each channel at a given redshift. Figure [3] 
shows the former for e + e _ pairs and photons, for the 
“best” approach. As a cross-check, we also display the 
results of the “approx” approach for e + e~ pairs using the 
“3 keV” baseline, re-adding the power lost to continuum 
photons for channel 5; the photon and “SSCK” cases are 
qualitatively similar. 

In general, we see that heating of the gas and distor¬ 
tions of the CMB spectrum dominate the energy losses 
at high redshifts prior to recombination, while excita¬ 
tion and ionization become more significant after recom¬ 
bination; in general, the fraction of power eventually de¬ 


posited falls for all channels at lower redshifts of injection 
(with some limited exceptions), as the universe becomes 
more transparent. Especially for e + e _ pairs, but also for 
photons, there is a striking structure at injection ener¬ 
gies around 10-100 MeV, where the fraction of injected 
power proceeding into ionization and excitation is quite 
high (and the fraction of power proceeding into heating 
is also enhanced); at lower energies, ~ 1 — 10 MeV, these 
channels are suppressed, and instead the production of 
low-energy continuum photons is enhanced. This accords 
with the behavior shown in Figure [2] As expected, the 
results of the “approx” and “best” methods are in general 
very similar for redshifts z < 1000. 


V. CORRECTED DEPOSITION-EFFICIENCY 
FUNCTIONS 

A. Deposition efficiency by channel 

Above we have presented results for the fraction of in¬ 
jected power deposited to the gas over the age of the 
universe; however, as discussed in Section |TTJ the key fig¬ 
ure of merit is instead the power deposited at any given 
redshift. Once an energy injection history is specified 
(as a function of redshift), one can use the results pre¬ 
sented above to integrate over the redshift of injection, 
and determine the deposited power originating from en¬ 
ergy injections at all earlier times. 

It is often convenient to normalize the deposited power 
at a given redshift to the injected power at the same red¬ 
shift (where both quantities are defined within a given 
comoving volume, per baryon, etc). This ratio defines 
an effective deposition efficiency curve f(z). Since de¬ 
position at a given redshift may have contributions from 
power injected at much earlier redshifts, f(z) can in some 
cases be greater than 1, but it is typically 0(0.1 — 1) [12]. 4 

As discussed above, it is not in general sufficient to 
derive the deposition-efficiency curve and then multiply 
the result by a model-independent prescription y[? ase (z) 
for the fraction of deposited power proceeding into the 
various deposition channels. Instead, we will define f c (z ) 
curves corresponding to the individual channels, which 
give the power deposited at a given redshift to a specific 
channel , normalized to the total injected power at the 
same redshift. By definition, f(z) = Ylcfc( z )- 

Given an injection of some species with redshift and 
energy dependence such that the rate of particle injec¬ 
tion per unit time per unit volume is given by d ^^ dt = 


4 In the context of annihilating DM, f(z) curves are usually de¬ 
fined with respect to the injected power from the smooth DM 
distribution, since this is easy to characterize, even though the 
onset of structure formation may greatly increase the injected 
and hence deposited power. In this paper we will not take these 
effects into account, and so will always define /(z) curves with 
respect to the total injected power. 
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FIG. 3: Total power deposited (down to 1 + 2 = 10) into each channel, as a function of injection energy and redshift. From top 
to bottom, the rows correspond to H ionization, He ionization, Lyman -a photons, heating, and sub-10.2 eV continuum photons. 
The left column describes energy absorption for e + e“ pairs (the a>axis “energy” label here indicates the kinetic energy of a 
single member of the pair at injection), while the right column describes energy absorption for photons. The center column is 
an alternate method of estimating the same quantities shown in the left column, as described in Section [IV A 1| and should be 
regarded as a cross-check. 
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I(z,E ), the corresponding f c {z) curves can be approxi¬ 
mated by: 

, _ Ej Efc E j I(z k , Ei)dV{z k )dt(z k )T^ es dE j 
}c ^ Z ’ ~ J2j EH{z\ Ei)dEHV{z i )dt{z i ) 

H{z i ){l + z*) 3 
~ V, E')dE’ 

X V 7 - * , , V E j I(z k ,&)T sp T s dEi. 

^ (1 + z k ) 3 H(z k )^ K 1 c ’ ljk 

k K 7 3 

( 8 ) 


Here the indices i, j, k label redshift of deposition, energy 
of injection and redshift of injection, as above; dt(z ) is 
the time interval corresponding to dln(l + z), and we 
have employed the relations H{z) = — dln(l + z)/dt and 
dV (zi) / dV{z 2 ) = (1 + £ 2) 3 / (1 + 2 i) 3 - 

For example, for conventional DM annihilation, with 
a rate that scales as the square of the density, I(z,E) oc 
(1 + z) e n(E), with the remaining factors being redslrift- 
and energy-independent. For DM decay with a lifetime 
much longer than the age of the universe, I(z,E) oc 
(1 + z) 3 n(E). In both cases n(E) describes the spectrum 
of injected particles for the species in question. Substi¬ 
tuting these expressions into Equation [8] and summing 
over all deposition channels reproduces the results of |I2| 
(up to numerical error associated with the discretization 
of the transfer function T). If we consider only particles 
injected at a single energy, so I(z,E) is proportional to 
a delta function in energy at , the / c (~) curves can be 
simplified to: 




h(z>)( 1 + ^) 3 v - Ej K%t es ,. 

I{z\Ei) ^ (1 + z k ) 3 H(z k )' U 


More generally, for any energy injection history I(z,E) 
that is a separable function of z and E, i.e. I(z,E) = 
I(z)n(E), the f c (z) curve can be written in the form: 


fctf) = 




( 10 ) 


Thus characterizing the f c (z ) curves for individual ener¬ 
gies is sufficient to describe all separable energy injection 
histories. We plot these curves for a range of injection 
energies, for the injection profile corresponding to DM 
annihilation, in Figure |4j (As a default, we present re¬ 
sults based on the “best” method described earlier, but 
as a cross-check, we also show results for the “approx” 
method for injected e + e _ and the “3 keV” baseline pre¬ 
scription.) 

By replacing T®^ les in Equation |s] with j^®= les ; base 
as defined in Equation [7J we can construct alternative 
curves, which correspond to rescaling the over¬ 
all f(z) by the channel-dependent y)? ase (z) factors and a 
correction factor to account for the losses into continuum 
photons. From Equation [7] it follows that, 

fr(z) = [f{z) f COII (z)} x h c ase (z), (11) 


where f CO rr(z ) is obtained by replacing T apa< £ es ’ base with 
in Equation g 

We can now define new, model-dependent fractions 
Xc(z) by Xc(z) = fc{z)/f(z)-, i.e. the fraction of de¬ 
posited power proceeding into each of the deposition 
channels. Likewise we define Xcorr(z) = fcon(z)/ f(z). 


B. Correcting the f(z ) curves for use with earlier 
studies 

Since the f c {z) are the quantities that determine the 
observable effects of energy injection (ionization, heat¬ 
ing, etc), a correct computation of the constraints us¬ 
ing the f(z) deposition-efficiency curves would employ 
the Xc{z ) fractions derived here. However, many con¬ 
straints have already been set assuming older, model- 
independent forms for the Xc ase (z) fractions. Since only 
the product f{z)Xc{z) matters, the use of incorrect Xc(z) 
fractions can be compensated by a correction to the 
f(z ) deposition-efficiency curve (this observation was 
also made in [Ill 120] ). However, since f(z) is channel- 
independent, only the power deposited to a single channel 
c can be completely described in this way; the deposition 
to other channels will be only approximate. Fortunately, 
this is not a problem for constraints that depend almost 
entirely on a single deposition channel. For example, 
constraints from the CMB anisotropies are primarily de¬ 
termined by ionization; while in principle the H and He 
ionization contributions should be treated separately (as 
in e.g. 0), in practice the contribution from He ioniza¬ 
tion is negligible, and so we can simply consider the sum 
of channels 1 and 2. 

When the key figure of merit for a particular constraint 
is set by f c (z) for some channel c, in order to adapt 
older studies performed using some “base” prescription 
for the Xc ase (£) fractions, one should define a new cor¬ 
rected deposition-efficiency curve J c ’ base (~) by: 

f cMse {z) = f( z ) Xc ( Z ) /x h c aSB ( Z ), (12) 

so that Xc ase (^)/ C,base (^) = fc{z). 

The alternate simplified method described above for 
computing the energy deposition fractions by channel, 
i.e. rescaling a simplified “base” prescription to ac¬ 
count for losses to continuum photons, corresponds to us¬ 
ing a channel-independent corrected deposition-efficiency 
curve f slm (z) given by: 



= f{z) - fcorr(z) 

= f 0) [1 - Xcorr(z)] • (13) 

Note that in this prescription the correction factor to 
f(z) is independent of channel (for channels 1-4), and is 
also independent of the choice of “base” prescription to 
describe the energy deposition by channel. This means 
that slightly different constraints will be obtained if this 
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FIG. 4: Power absorbed into each channel from particles injected by DM annihilation (or another process scaling as cosmological 
density squared), as a function of injection energy and redshift of absorption, normalized to the total injected power at the same 
redshift. From top to bottom, the rows correspond to H ionization, He ionization, Lyman-a photons, heating, and sub-10.2 
eV continuum photons. The left column describes energy absorption for e + e~ pairs (the x-axis “energy” label here indicates 
the kinetic energy of a single member of the pair at injection), while the right column describes energy absorption for photons. 
The center column is an alternate method of estimating the same quantities shown in the left column, as described in Section 
|IV A 1| and should be regarded as a cross-check. 
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FIG. 5: Corrected f(z) curve for particles injected by DM annihilation, as a function of injection energy and redshift of 
absorption. In the left panel we use the “3 keV” baseline ionization fractions (so these f(z) curves should be used with analyses 
that employed the same prescription); in the center panel we use the “SSCK” baseline. In the right panel we correct for the 
continuum losses using the results of Figure|hj and thus derive an alternate channel-independent f Blln (z) curve. The upper row 
describes e + e“ pairs (the ir-axis “energy” label here indicates the kinetic energy of a single member of the pair at injection), 
the lower row describes photons. 


simplified prescription is combined with studies using 
different “base” prescriptions; in contrast, if the cor¬ 
rect J c . base ( 2 ) curve is employed, the dependence on the 
“base” prescription will cancel out between the corrected 
f(z) curve and the choice of the Xc ase ( 2 ) factors in the 
original analysis. 

In Figure [ 5 ] we plot the j lon . base (~) curves for the 
“SSCK” and “3 keV” choices of “base” prescription, 
again for an annihilation-like history: these curves con¬ 
stitute our best estimate of the appropriately corrected 
deposition-efficiency curves for the purposes of comput¬ 
ing CMB constraints on DM annihilation. We also dis¬ 
play the / Slm (z) curves obtained as described in Equa¬ 
tion |13| In Figure [6] we plot the approximate correc¬ 
tion factor Xcorr (z), which should be interpreted as the 
fraction of deposited energy proceeding into previously 
unaccounted-for continuum photons, for the energy in¬ 
jection history corresponding to conventional DM anni¬ 
hilation. 

From Figure[6]we see that the correction to f(z) due to 
continuum losses is largest at injection energies around 
1 — 100 MeV (depending on redshift) for photons, and 
at slightly lower energies (~ 1 — 10 MeV) for e + e - 
pairs. This is consistent with the discussion of Figure 
[5] the correction is smaller than one might expect for 
non-relativistic e + e _ pairs (with injection kinetic ener¬ 
gies well below 1 MeV) because most of the injected en¬ 


ergy is bound up in their mass, and thus the deposition 
of the kinetic energy is almost irrelevant. Such particles 
deposit their energy primarily through the annihilation 
of the positron, producing photons at 511 keV and below. 

The correction factor falls abruptly at both lower and 
higher energies, although it is still appreciable at the 
highest energies. The correction factor for high-energy 
injected particles is somewhat complex, as it depends 
on the discontinuous cooling of the high-energy photons 
and electrons. For example, a 5 GeV electron injected 
at 2 ~ 1000 will upscatter a ~ 1 eV CMB photon to an 
energy of order q 2 ~ 10 8 eV; the resulting photon will 
dominantly lose its energy by pair production or Comp¬ 
ton scattering [T2j, partitioning its energy into lower- 
energy electrons and positrons. The upper end of this 
spectrum, corresponding to 10-100 MeV electrons, will 
again upscatter CMB photons, now to energies ranging 
from hundreds of eV to tens of keV; these photons will 
efficiently ionize the gas. (All numbers in this paragraph 
are approximate and for illustration only.) An order of 
magnitude reduction in initial electron energy would re¬ 
duce the typical energy of upscattered photons by two 
orders of magnitude, to around 1 MeV; the resulting sec¬ 
ondary electrons would upscatter CMB photons to ener¬ 
gies too low to further interact with the gas, thus losing 
a large fraction of their energy to the continuum. How¬ 
ever, another order of magnitude reduction in the initial 
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FIG. 6: Fraction of deposited energy proceeding into previously-unaccounted-for continuum photons, y CO rr( 2 ), for a DM- 
annihilation-like injection history. The upper panel gives results for e + e“ pairs (the rr-axis “energy” label here indicates the 
kinetic energy of a single member of the pair at injection), the lower panel for photons. 


electron energy would mean that CMB photons would 
be upscattered to 0(10 keV) energies, at which point (at 
z = 1000) they would be efficient photoionizers them¬ 
selves and produce electrons in an energy range where 
atomic processes dominate the cooling. Patterns of this 
type are the reason for the “striping” visible in Figures 

m 


Specifically, we define these curves by: 

t g(^)(l + ^) 3 v 1 

/i°ssG J ^ E3l(z\ Ei)dEi ^ (1 + z k ) 3 H{z k ) 

x £ WI{z\&)SZ s %dE\ (14) 


f H(z l )(l + z 1 ) 3 _1_ 

thigh, cG ) Eil{z\ Ei)dEi ^ (1 + z k ) 3 H(z k ) 

x Y, Eil ( zk > Ej ) S l%t S dE j , (15) 

C. Low-energy particle spectra from arbitrary :l 

energy injections 


In the same spirit as the f(z) curves, which describe 
total absorbed power, one can integrate over injection 
redshift to determine the low-energy electron and pho¬ 
ton spectra produced by a specific energy injection his¬ 
tory, at any redshift. Likewise, one can produce a fioss(z) 
curve describing depletion of the CMB spectrum by scat¬ 
tering, and /high,c (z) curves that include only the power 
deposited by the cooling of high-energy electrons, from 
the initial outputs of the code described in Section m 
Such results can be converted into the f c (z) curves, us¬ 
ing results from codes that describe the cooling of low- 


energy electrons, just as described earlier in Section IV 


the photon spectrum at very low energies may also con¬ 
stitute an observable in its own right, as a distortion to 
the CMB blackbody spectrum. For convenience we can 
normalize the photon and electron spectra to the number 
of injected pairs at the “output” redshift (when the low- 
energy electrons and photons are produced), thus can¬ 
celing out model-dependent normalization factors in the 
energy injection rate; we denote these normalized low- 
energy spectra by F sec (z, E sec ). (Note this choice is not 
identical to the f[z) curves, where we normalize to in¬ 
jected power rather than number of annihilations.) 


F sec (z\E l se c ) 


fJ(2 i )(l + ^) 3 1 

'E j I(* i ,E3)dE3 “ (1 + z k ) 3 H(z k ) 

(16) 

3 


For convenience, we include these arrays in our pro¬ 
cessed results, for the energy injection history corre¬ 
sponding to conventional DM annihilation (i.e. a fixed 
spectrum of annihilation products and a I(z,E) oc (1 + 
z) 6 redshift dependence). 


VI. CONCLUSIONS 

We have presented the results of a comprehensive nu¬ 
merical study of the energy losses of keV-TeV photons 
and e + e~ pairs in the cosmic dark ages, from their in¬ 
jection energies down to the ~ 3 keV scale where inter¬ 
actions with the gas begin to dominate the energy losses 
of electrons. All our results are provided in .fits and 
.dat format to facilitate the matching of these results 
onto detailed models of the low-energy cooling. We have 
employed previously published results for the low-energy 
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cooling to estimate the partition of deposited energy be¬ 
tween ionization, heating and the production of Lyman- 
a and continuum photons. We have demonstrated how 
to use these results to compute corrected deposition- 
efficiency f(z) curves for use with studies of constraints 
on energy injection, and similarly made these processed 
results public. 
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Appendix A: Summary of output files 

We provide three types of output files: resultsgrid_species_fulloutput.fits (where “species” can be either 
“phot” or “elec”, corresponding to injection of photons and e + e~ pairs respectively), species.results.derived. fits 
(likewise), and supplementary_deposition_fractions.fits. The first describes the direct outputs of our code - 
energy deposited by electrons above 3 keV, and low-energy electron and photon spectra produced in each timestep 
whereas the second describes the processed results, giving our estimate of the energy deposited into each of the 
five channels described in Section |III[ The third file provides the reference ionization history we used, the resulting 
“SSCK” fractions, and the table of \c fractions as a function of redshift and electron energy derived from [[!]■ The 
contents of each type of file are as follows: 


species.data.fits 


• OUTPUT_REDSHIFT: this 63-element array provides the abscissa for output redshift, i.e. the value of l+z at which 
the energy is deposited. 

• ENERGY: this 40-element array provides the abscissa in energy, with values given by log 10 (energy in eV). Note 
that this is kinetic energy of one of the two particles in the case of e + e _ pairs; a particle annihilating or decaying 
to e + e - would need a mass sufficient to provide this energy in addition to the mass energy of the pair. 

• INPUT .REDSHIFT: this 63-element array provides the abscissa for input redshift, i.e. the value of 1 + z at which 
the energy is injected. 


• CHANNELS : this 5-element array lists the five deposition channels: hydrogen ionization, helium ionization, Lyman- 
ol / excitation, heating and continuum photons. 

• DEP0SITI0N_FRACTI0NS_0RIG: this 63x40x 63 array provides the table T^ ecies = Yl, c f° r appropriate 

species: that is, for a particle injected at some input redshift and energy (given by the abscissa arrays), the 
fraction of its initial energy deposited to all channels in the (log-spaced) timestep associated with the output 
redshift. Up to small numerical differences, this table should match the DEPOSITION_FRACTIONS table in the 
files associated with m 

• F.ORIG: this 63 x 40 array describes the original depositon-efficiency /(t)-curve for DM annihilation to the species 
in question, with DM mass given by ENERGY, sampled at the redshift points given by the OUTPUT .REDSHIFT array. 
This is based on the total energy deposition and contains no corrections to account for model-dependent Xc{z) 
fractions. 


• CMBL0SS.FRACTI0NS: this 63 x 40 x 63 array, labeled 5^®“®*. in the text, describes the power scattered out of the 
CMB, as a fraction of the initial energy of the injected particle, in the timestep corresponding to the deposition 
redshift. Its index structure is the same as that of DEP0SITI0N.FRACTI0NS.0RIG. 


HIGHDEP.FRACTIONS: this 63 x 40 x 63 x 5 array provides the table <S'® p 7^ les , as defined in Table |l| that is, the 
power deposited into each of the channels c = 1 — 5 by cooling of high-energy electrons (see Section III ), as a 
fraction of the initial energy of the injected particle, in the timestep corresponding to the deposition redshift. 
This array’s first three indices are the same as those for DEP0SITI0N_FRACTI0NS.0RIG; the last index corresponds 
to c, and lists the 5 different channels. As noted in Section ITTT 


channel 2 is always empty. 
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• PHOTENG: this 40 x 500 array provides the abscissa for the low-energy photon spectrum, the energy in eV for 
the low-energy secondary photons (500 energy bins), at each of the 40 injection energies defined by the ENERGY 
array. 

• ELECENG: this 40 x 500 array provides the abscissa for the low-energy photon spectrum, the kinetic energy in 
eV for the low-energy secondary electrons (500 energy bins), at each of the 40 injection energies defined by the 
ENERGY array. 

• L0WENGPH0T_SPEC: this 63 x 40 x 63 x 500 array provides the table as defined in Table |lj that is, the 

spectrum dN/dE of low-energy photons produced per injected pair of particles, in the timestep corresponding 
to the deposition redshift. This array’s first three indices are the same as those for DEP0SITI0N.FRACTI0NS; the 
last index l corresponds to the energy of the secondary photons, defined by the (j,l )th element of the PHOTENG 
array. 

• L0WENGELEC_SPEC: this 63 x 40 x 63 x 500 array provides the table S h p p ^// S . as defined in Table |lj that is, the 
spectrum dN/dE of low-energy electrons produced per injected pair of particles, in the timestep corresponding 
to the deposition redshift. This array’s first three indices are the same as those for DEP0SITI0N_FRACTI0NS; the 
last index l corresponds to the kinetic energy of the secondary electrons, defined by the (j, Z)t.li element of the 
ELECENG array. 


species_processed_results.fits 

• OUTPUT .REDSHIFT: this 63-element array provides the abscissa for output redshift, i.e. the value of l + z at which 
the energy is deposited. 

• ENERGY: this 40-element array provides the abscissa in energy, with values given by log 10 (energy in eV). Note 
that this is kinetic energy of one of the two particles in the case of e + e _ pairs; a particle annihilating or decaying 
to e + e~ would need a mass sufficient to provide this energy in addition to the mass energy of the pair. 

• INPUT .REDSHIFT: this 63-element array provides the abscissa for input redshift, i.e. the value of 1 + z at which 
the energy is injected. 

• CHANNELS: this 5-element array lists the five deposition channels: hydrogen ionization, helium ionization, Lyman- 
a, heating and continuum photons. 

• DEPOSITIQN.FRACTIQNS.ORIG: as above. 


• F.ORIG: as above. 

• DEP0SITIQN.FRACTI0NS.NEW: this 63 x 40 x 63 x 5 array provides the T®^ les array, summing the contributions 
from high-energy deposition and integrating over the low-energy photon and electron spectra. This array’s first 
three indices are the same as those for DEPQSITIQN.FRACTIONS.QRIG defined above; the last index corresponds 
to c, and lists the 5 different channels. 

• C0NT.CDRR.FRACTIDNS: this 63 x 40 x 63 array provides the T® oiiTjk arra y> the contribution to channel c = 5 
arising from the low-energy photons produced at each timestep by cooling of electrons above 3 keV (divided as 
usual by the initial injected energy, and having corrected for the original energy of these low-energy photons). 
This array’s indices are the same as those for DEP0SITI0N.FRACTI0NS.0RIG defined above. 

• F.I0N: this 63x40 array describes the deposition-efficiency /-curve for DM annihilation to the species in question, 
corrected to give the true power into ionization combined with the “3 keV” prescription. This quantity is denoted 
jion,3keV^ j n Section [v] The DM mass is given by the ENERGY array, and the function is sampled at the redshift 
points given by the 0UTPUT.REDSHIFT array. 

• F.C0RR: this 63 x 40 array describes the deposition-efficiency /-curve for DM annihilation to the species in 
question, approximately corrected by rescaling the deposited energy according to the unaccounted losses into 
continuum photons. This quantity is denoted f slnl (z) in Section [v] The DM mass is given by the ENERGY array, 
and the function is sampled at the redshift points given by the 0UTPUT.REDSHIFT array. 


F.CMBLDSS: this 63 x 40 array describes the power scattered out of the CMB at each redshift, as a fraction of 
the power injecte d at that redshift, for DM annihilation to the species in question. This quantity is denoted 
/ioss(~) hi Section V C The DM mass is given by the ENERGY array, and the function is sampled at the redshift 
points given by the 0UTPUT.REDSHIFT array. 
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F.HIGHDEP: this 63 x 40 x 5 array describes the power deposited to each channel by the cooling of high-energy 
electrons at each redshift, as a fraction of the power injected at that redshift, for DM annihilation to the species 
in question. This quantity is denoted /high,c ( 2 ) in Section V C The DM mass is given by the ENERGY array, and 
the function is sampled at the redshift points given by the OUTPUT .REDSHIFT array. 


• PHOTENG: this 40 x 500 array provides the abscissa for the low-energy photon spectrum, the energy in eV for 
the low-energy secondary photons (500 energy bins), at each of the 40 injection energies defined by the ENERGY 
array. 


• ELECENG: this 40 x 500 array provides the abscissa for the low-energy photon spectrum, the kinetic energy in 
eV for the low-energy secondary electrons (500 energy bins), at each of the 40 injection energies defined by the 
ENERGY array. 


FSPEC.PH0T: this 63 x 40 x 500 array provides the spectrum cLN/dE of low-energy photons produced at each 
redshift, normalized to the number of DM annihilations occurring at that redshift, for DM annihilation to the 
species in question. This quantity is denoted F' y (z, E-y) in Section V C The DM mass is given by the ENERGY 
array, and the function is sampled at the redshift points given by the 0UTPUT.REDSHIFT array; the energies of 
the secondary photons are given by the corresponding elements of the PHOTENG array. 


• FSPEC.ELEC: this 63 x 40 x 500 array provides the spectrum dN/dE of low-energy electrons produced at each 
redshift, normalized to the number of DM annihilations occurring at that redshift, for DM annihilation to the 
species in question. This quantity is denoted F e (z,E e ) in Section VC The DM mass is given by the ENERGY 
array, and the function is sampled at the redshift points given by the OUTPUT .REDSHIFT array; the kinetic energies 
of the secondary electrons are given by the corresponding elements of the ELECENG array. 


deposition.fractions.supplement.fits 

• REDSHIFT: this 63-element array provides the abscissa for deposition redshift. 

• XH: this 63-element array provides the hydrogen gas ionization fraction as a function of redshift, for the baseline 
ionization history with no energy injection (calculated using RECFAST, as in [34]). 

• CHANNELS: this 5-element array lists the five deposition channels: hydrogen ionization, helium ionization, Lyman- 
a , heating and continuum photons. 

• SSCK: this 63 x 5 array describes the fraction of deposited power proceeding into channels 1 — 5 under the simple 
“SSCK” prescription, as a function of redshift (channels 2 and 5 are not populated by this prescription), i.e. 

xf CK (*)- 

• ELECTR0N.ENERGY: this 6-element array provides the kinetic energy values at which the CHANNEL.FRACTI0NS 
array is evaluated. 

• CHANNEL.FRACTI0NS: this 6 x 63 x 5 array describes the fraction of deposited power proceeding into channels 
1 — 5 in a detailed calculation of the low-energy cooling (presented in IJ3] and based on H3), for the six electron 
injection energies listed in ELECTR0N.ENERGY, as a function of redshift. 

The files are available online in .fits format, at http://nebel.rc.fas.harvard.edu/epsilon; we supply a 
Mathematica notebook demonstrating how to read the .fits files and reproduce the calculations in this note. Finally, 
we also provide the key results as . dat files, as described below: 

Overall best estimate for the energy deposition fractions: 

• species_channel=c_deposition_fractions.dat: holds T^®^ les for each choice of c and species (“phot” cor¬ 
responds to photons, “elec” to e + e~ pairs). Starting with the third column, the entries in the first row give 
the injection redshift 1 + z. Starting with the second row, the first column lists the deposition redshift 1 + z, 
and the second the injection energy (as defined in the ENERGY array in the .fits files); subsequent columns list 
T c ,ijk for the appropriate triple of injection redshift, deposition redshift and injection energy. (Holds the same 
information as DEP0SITI0N_FRACTI0NS_NEW in the .fits files.) 

• species.continuum.correction.fractions.dat: as species_channel=cJiighE.deposition.fractions.dat, 
except T^T replaced with T^^ k . (Holds the same information as C0NT.C0RR.FRACTI0NS in the .fits 
files.) 
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Raw outputs of the high-energy code: 

• species_channel=c JiighE_deposition_f ractions . dat: as species_channel=cTiighE_deposition_fractions . dat, 
except T-— is replaced with S^ p ^ es . (Holds the same information as HIGHDEP.FRACTIONS in the .fits files.) 

• species_cmbloss_fractions.dat: as species_channel=c JiighE_deposition_fractions. dat, except 

is replaced with (Holds the same information as CMBLOSS_FRACTIONS in the .fits files.) 

• species_lowEsecspectra.dat: holds for each choice of “sec” (“electron” for secondary electrons and 

“photon” for secondary photons) and each injected species. Starting with the fourth column, the entries in the 
first row give the injection redshift 1 + z. Starting with the second row, the first column lists the deposition 
redshift 1 + z, and the second the injection energy (as defined in the ENERGY array in the .fits files), the third 
the secondary (kinetic) energy, expressed as logio(-E se c/eV); subsequent columns list S f° r appropriate 
quadruple of injection redshift, deposition redshift, injection energy and secondary energy. (Holds the same 
information as L0WENGELECJ3PEC and LOWENGPHOT.SPEC in the .fits files.) 

f(z) curves for DM annihilation: 

• species_bestf_DMann_3keV.dat: holds j lon > 3keV ( 2 r), our best-estimate corrected f(z) curve for DM- 
annihilation-like energy-injection histories, for studies which assume the “3 keV” baseline prescription to set 
constraints on energy injection via the ionization channel (e.g. the recent constraints on DM annihilation pre¬ 
sented by the Planck Collaboration [33] 1. Starting with the second column, the entries in the first row give the 
deposition redshift 1 + z. Starting with the second row, the first column lists the injection energy (as defined in 
the ENERGY array in the .fits files); subsequent columns list f(z) for the appropriate pair of injection energy 
and redshift. (Holds the same information as F_I0N in the .fits files.) 

• species.simplef _DMann. dat: holds / slm (z), our simplified approximate prescription for the corrected f{z) curve 
for DM-annihilation-like energy-injection histories. Layout is the same as species_bestf_DMann_3keV.dat. 

(Holds the same information as F.CORR in the .fits files.) 
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